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INTRODUCTION 


The nature of the stress field in front of a crack lying in 
the surface between bonded dissimilar materials was first investi- 
gated by Williams (ref. 1) . He observed that not only do the 
stresses grow at a rate inversely proportional to the square root 
of distance from the crack tip, they also exhibit an oscillatory 
singularity with wave length inversely proportional to the absolute 
value of the logrithm of distance from the crack tip. The problem 
of calculating stress intensity factors for various special load- 
ings and geometries has been treated by other authors, among them 
Erdogan (ref. 2 and 3), England (ref. 4), Rice and Sih (ref. 5), 
and Erdogan and Gupta (ref. 6) . In all cases for which results 
are given the region is unbounded and the loads are uniform. 

For more general boundary value problems involving imperfect 
bonding of dissimilar materials numerical methods must be resorted 
to, and both the growth and oscillatory nature of the singularity 
can be expected to cause numerical difficulties. In addition, be- 
cause the elastic moduli of the materials are generally different, 
discontinuities in components of stress and strain develop nat- 
urally on the bond. Recently we extended the contour integral 
method to problems of this type. It turns out that the nature of 
the loading and restraints , even on remote edges , can have a sig- 
nificant effect on the stress intensity. In this paper we treat 
some example problems to illustrate this. 


This work was supported in part by a grant from the National 
Science Foundation. 
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CONTOUR INTEGRAL METHOD 


The basic boundary value problem is illustrated in Fig. 1. 

Two dissimilar materials are joined along a straight edge with one 
or more cracks present. The composite is loaded or restrained on 
the remote boundary, and the crack faces are free of load. Local 
cartesian and polar coordinates are introduced with origin at a 
crack tip and the negative x-axis (0 = ± it) along the crack edges. 
The subscript 1 is arbitrarily assigned to material below the axis 
(— 7T £ 0 £ 0) , and the subscript 2 is used for the other material 
(0 < 0 < it) . Also introduced is the so-called bimaterial constant 


^1 + ** 2*1 
y 2 + u 1 k 2 


( 1 ) 


where y^, y ^ are the respective shear moduli and k = 3-4v for plane 
strain, or k = (3-4v) (1+v) for plane stress, v being Poisson's Ratic 

Notation for complex displacement and stress fields in terms 
of components referred to the local polar coordinate system are 
introduced as follows: 


u = u r + iu e 


°r = °rr + iT rS 


a e a ee 


- IT 


(2) 
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Then the displacement and stress fields in the neighborhood of the 
crack tip in each material are of the form-2 


u, = 


K r X [K 1 * 1 ? X “ 1)0 - ye i( - X - 1)e ] 


1 " 2y 1 (1+y) X * L "l 


+ . ■ : *r~ r X [e l( X ^ 4 * * * * 9 - e 1 ^ X+1)9 ] + remainder 


2y 1 (l+y) 


u - j r X [ K Te i(X ~ 1)9 + e i( " X-1)e i 

U 2 “ 2y 2 (l+y)X r lK 2 Y ® + e J 

+ r X fe i( - X+1)9 - e i( - X+1 > 0 l + remainder 


lr 


k x-i, i(x-i)e _ iH-De, 

(1+y) L Y J 


2r 


+ jjjyf r X-1 (Xe l(X 1)9 - (X-2)e l( X+1)9 ] + remainder 


K r X-l f i(X-l)0 _ i(-X-l) 0, 

(1+y) r lye e J 


(3) 


, Ky X-l r T i(-X-l)0 . i(-X+l)0, . . , 

+ ^ x+y ' ) r CXe - (X-2)e ] + remaxnder 


10 


4 K .... r X-l f i(X-l)0 i(-X-l)0, 

(1+y) r Le ye J 


+ r X “ 1 (e l( X+1)9 - e 1 (" X “D e ] + remainder 


’20 “ TI+yT 


K r X - 1 [ye i(X " 1)9 + e i(_X “ 1)9 ] 


+ -g&T r X - 1 [e i( ' X+1)9 - + remainder 


'Except for notational differences these results were also obtained 
in references 3, 4 and 5. 
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where 


A = | + I? ln V = I + ie (e = h ln Y) (4) 

i B 

and K = K Q e p is a complex stress intensity factor with the 
following "physical" interpretation: 


lim a 
r-+0 


9 


9=0 


= -lim (a yy - iT xy ) 
x->0 


y=0 


= Kr A_1 = K r 'V (eln r + 6) 
o 


hence on the bond immediately in front of the crack tip we have 


, - 5>. 

yy /f- 

K 


xy 


cos (eln r + 3) + remainder 


sin (e In r + 3) + remainder 


(5) 


/r 


Thus K q governs the amplitude growth rate of both the normal stress 

and shear stress while 3 determines a nonsignificant . phase shift. 
The complex crack opening displacement is also governed by the 


COSh £7T 7 } A 

TT+yPT k r 


stress intensity factor: 



Au = u 0 


* 

«1 + 


0 = 7T 

0=~7T 

^2 

v * 


The amplitude of the complex crack opening displacement can be put 
in the form 


Au 


< Ul + VjKj) (Vj + Vl * 2 ) cosh ET 
yiP 2 (Vl + u 2 + v 2 < 2 + /t — — 2 

/ 4 + e 


K/r~ 

O 


( 6 ) 


A contour, integral representation for the stress intensity 
factor is obtained from the reciprocal work identity by introduc- 
ing a suitable artificial singular elastic state. Briefly, 
we observe that for arbitrary values of the complex constant C, 
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a singular elastic state corresponding to zero body force and with 
no traction on the lines 0 = ± it is defined in the bimaterial 
region by 


2y ] _u* 


2U 2 »* 





CXr" X [e i(X+1)0 - e 1(X ‘ 1)9 ] 

+ Cr- X [ Kl e 1( - X - 1)e - 

C'vr" ' [e 1 * ' +1 ' 0 - 
+ cr-^l^ye 1 '-^ 119 - e i<X - 1)0 ] 

CAr 1 IXe 1 ^ 9 - (A+2)e 1 * X+11 6 ] 

+ CX J :- X - 1 [Ye i<X - 1)e - e i( - X - 1)0 ] 

C Y Ar~ X " 1 [Ae :L(X_1)e - U+2) e 1 (X+1) 6 ] 
+ CAr-^he 1 '*- 1 ’ 6 - Ye 1 '" 1 ' 1 ’ 6 ] 

aA 2 r- X - 1 [e i(X+1 > e - e i(X - 1)6 ] 

- CAr' X ~ 1 (e i( ' X_1)0 + Y e ilX ~ 1>9 I 


(7) 


a e ; = CYA 2 r - X - 1 te 1(X+1 > 9 - e 1 *^ 1 * 9 ] 

- CAr-^lYe 1 '- 1 - 116 + e i(X - 1)e ] 

This elastic state has the further property that on the contour 
C £ (a circle of radius e centered on the origin) we calculate a 

finite contribution from the reciprocal work as the contour 
shrinks to a point: 
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lim (u**t - u*t ) ds 

e - o 4 ~ ~ 


im 

>0 J 


lim 

e-*0 


— * — * 

Re [u cr - ua ] ds 
r r 


-ir(y 1 + y^) 


Re CK 


( 8 ) 


Upon noting that the reciprocal work vanishes on the complete 
contour C Uc + Uc Uc indicated in Pig. 1 as a consequence of Betti's 

° ^ -f- ' — 

theorem, and on the crack edges C UC since the tractions in any 
case vanish there, we obtain the representation 


I Re Kl 
\lm K I 



* 

u 


• t) ds 


(9) 


y x y 2 

where for Re K we choose C - - — r- — y in Eq. (7) in calcu- 

~ + ^ 2 k 1 ' 11 11 

* * x 4 x y,y 2 

lating t and u , whereas for Im K we take C = -i — r- . 

~ ~ 7T 1 + ^2 k 1 

The values of u and t on the contour C Q are obtained numerically. 

For the results given in this paper we used code TEXGAP (ref. 7) 
which performs isotropic linearly elastic plane analyses using 
conventional quadratic displacement triangles and isoparametric 
quadrilaterals . 


NUMERICAL RESULTS 

The four cases treated involve a finite bimaterial strip 
loaded in tension and are sketched in Fig. 2. From symmetry con- 
siderations we need to consider only the shaded region, and in 
Fig. 3 we show a typical grid (symmetrically defined in each half 
region) for the finite element analyses. Half the contour used 
for evaluation of the stress intensity factors is shown in dashed 
line in Fig. 3. We note that the four distinct problems con- 
sidered are obtained from the same grid and boundary conditions 
on the edges parallel to the crack, but with the following boun- 
dary conditions on the vertical edges: 
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i) Central crack - free edges: AB restrained, CD un- 

restrained 

ii) Central crack - fixed edges: AB and CD restrained 

iii) Single edge crack: AB and CD unrestrained 

iv) Double edge crack: AB unrestrained, CD restrained. 

Two sets of results are plotted in Pig. 4 and 5. The first 
shows the effect of different crack sizes in a given strip for 
each case; the second shows' the effect of changing the relative 
dimensions of the strip for a fixed crack length to width ratio. 

In each case the results are normalized using the stress intensity 
factor for an infinite region loaded uniformly in tension normal 
to the crack and restrained from motion parallel to the crack on 
the remote boundary. This case is equivalent to an infinite bi- 
material plate with vanishing stresses at infinity and a uniformly 
pressurized crack on the bond, for which analytical results are 
given in references (4) and (5) : 

K = A + 4e a_/si/2 (10) 

00 Q 

It is interesting to note that for real materials the bimaterial 
constant y is restricted to values between 1 and 3; consequently, 
the maximum variation in that can be achieved by varying the 

properties of the two materials (this enters only through the 
parameter e) is less than six percent, thus the isotropic case 
furnishes an excellent (lower bound) estimate for K^. The data 

plotted in Fig. 4 and 5 are based on material properties 

K 1 = 1 . 6 , <2 = 1 * 8 / ^ 2^1 ~ ^ 

which yields the value y = 1.5. 
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